Revisiting the dynamical exponent equality z = d for the dirty boson problem 
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It is shown that previous arguments leading to the equahty z — d {d being the spatial dimen- 
sionality) for the dynamical exponent describing the Bose glass to superfluid transition may break 
down, as apparently seen in recent simulations (Ref. 8|). The key observation is that the major 
contribution to the compressibility, which remains finite through the transition and was predicted 
to scale as K ~ (where 5 is the deviation from criticality and u is the correlation length 

exponent) comes from the analytic part, not the singular part of the free energy, and therefore is 
not restricted by any conventional scaling hypothesis. 
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Beginning with the realization that the zero tempera- 
ture onset transition of superfluidity in a random medium 
should properly be treated as a fluctuation driven quan- 
tum phase transition P, [3, H, ^ number of scaling 
arguments have been put forward 0, Q to determine or 
to place bounds on the critical exponents describing it. 
Probably the most significant of these was an argument 
that the dynamical exponent, z, describing the relative 
divergence of the temporal and spatial correlation lengths 
via S,T should be equal to the dimension, d, of space. 

The argument was supported by exact calculations in 
d = 1, where indeed z = 1 0], byl-|-e renormalization 
group calculations Q which, however, are not completely 
rigorous due to the absence of a form for the Hamiltonian 
for noninteger d > 1, and by a series of quantum Monte 
Carlo (QMC) studies in d < 2 0, 0, H- Results in d > 2 
are restricted to a badly controlled double e-expansion 
which is inappropriate for testing an exact scaling argu- 
ment of this type 

m 

, but nevertheless points to values 
of z significantly larger than unity Q . 

With one controversial exception Q , the earlier QMC 
results were fit to z = d, but the numerical error bars 
were not very tight. The more recent QMC study in d = 
2 Q provided a much more stringent test: using larger 
systems and a joint, optimal critical functional form fit 
to three separate thermodynamic quantities, the value 
z = 1.40 ± 0.02 was found, strongly violating z = d. In 
light of this, the purpose of this paper is to revisit the 
scaling arguments [3, [3| that led to z = d, and show that 
they can in fact break down: z might be independent 
of all the other exponents. In the absence of further 
theoretical justification, we raise the possibility that z is 
unconstrained by any simple scaling argument. 

Before presenting details, we briefly review the essence 
of the scaling arguments and indicate where they could go 
wrong. The helicity modulus, T [or superfluid density, 
Ps = (m^/?i^)T], quantifies the response of the super- 
fluid to gradients in the phase, 0, of the order parame- 
ter. Specifically, the free energy density contains a cor- 
rection A/^ = {T/2PV) dr J d'^x|V(/)(x, r)p, where V 



is the volume, f3 — l/ksT, and r is the usual imagi- 
nary time variable. The essence of the Josephson scal- 
ing argument is that since Vc/) has dimensions of in- 
verse length, it should scale as ^ \S\'^ where 5 is 
the deviation of the control parameter from criticality 
and v is the correlation length exponent. On the other 
hand the singular part of the free energy, fs, is de- 
fined to scale as |(5p~". This implies that T ~ \S\'" 
with V ~ 2 — a — 2v = {d -\- z — 2)v, where the last 
equality follows from the quantum hyperscaling relation 
2 — a = (d-f z)v [3]. Now, the Josephson relation, which 
connects changes in the chemical potential to the time 
derivative of the phase, actually allows one to interpret 
the compressibility, helicity modulus in the imag- 

inary time direction. Thus k enters a free energy correc- 
tion A/^ = (At/2/31/) dr J d'^x{dr(l)Y. The same logic 
thenn suggests that dr cj) should scale as ~ |<5|^'^, lead- 
ing to K '--^ |(5|'"^ with Vt = 2 — a — 2zv = (d — z^v. Since 
both the superfluid and the Bose glass phases have finite 
compressibility, one expects on physical grounds that n 
should be finite and nonzero at i5 = as well. This im- 
mediately leads to the prediction z = d 0, 0]. However, 
a hidden assumption in these arguments is that A/a; and 
A/t- can be included in the singular part of the free en- 
ergy, /g. We shall show that this is correct for A/a; but 
not for A/.r- In fact the main contribution to A/^ comes 
from the nonsingular analytic part of the free energy, 
/a, so that K is dominated by its analytic part which is 
trivially finite at the transition, and /^ yields only small 
corrections to this that vanish at the critical point. 

Rigorous definitions of the helicity moduli! involve 
comparing free energies with twisted and untwisted 
boundary conditions. Previously such boundary condi- 
tion dependence. A/, has always been included in jg. 
normally the strong dependence required for a finite he- 
licity modulus requires long or quasi-long range order 
that is present only in the superfluid phase. The fact that 
T = in the disordered phase then guarantees that it can 
arise only from singular terms in the free energy. In the 
present problem, however, k is nonzero in both phases: 
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the gapless excitation spectrum of the Bose glass phase 
leads to power law order in imaginary time (though not in 
space) [3| . The free energy of both phases therefore have 
strong temporal boundary condition dependence, and an 
analytic contribution is actually very natural. 

From our analysis there emerges the following general 
criterion for when A/ can be dominated by fa- If the 
twisted boundary condition breaks a fundamental sym- 
metry of the model, then it can be expected to generate 
relevant (in the renormalization group sense) terms in 
the Lagrangian. This leads to a new relevant, diverging 
scaling variable in the singular part of the free energy 
that generally will dominate all other contributions to 
A/. In this case the usual Josephson scaling relation 
will hold. On the other hand, if no additional symmetry 
is broken, no new relevant scaling variable results and 
the twist will lead only to small shifts in the parameters 
already present in the untwisted Lagrangian. The helic- 
ity modulus, which involves derivatives with respect to 
these shifts, will then be dominated by fa- In the present 
problem the twist couples to particle-hole symmetry 
which is always broken at the Bose glass to superfluid 
transition. 

For convenience we consider a continuum |7/>|'*-model 
functional integral representation of the partition func- 
tion, Z — J D-if: ex:p{—^B), with Lagrangian ^, ^] 



dr / d'^xl - j7/;*(x,T)V2^(x,r) 



if^*(x,T)[9, -Ai(x)]2^(x,T) 

r(x)|i^(x,T)|2 +M|V;(x,r)|' 



(1) 



where J « fi? /2m > is the boson hopping amplitude 
which favors spatial ferromagnetic order in the phase of 
the field ip; /i(x) is the (static) random external potential; 
and K « 1/2mo > 0, where uq is the soft core repulsion 
between nearby bosons, favors ferromagnetic order along 
the imaginary time dimension as well. This model is 
a continuum approximation to the Josephson junction 
array Lagrangian in which ^/'(x, r) — > exp[i(/)i(T)] where 
0i(T) is the Josephson phase at site «, and the 'ip*V'^ip 
term is the continuum limit of the Josephson coupling, 
Jij cos[(/>i(T) — (j)j{T)], between two nearby sites i and j. 
The rjV'P + ulV'l^ terms then represent the usual soft 
spin approximation to the the constraint [ipl = 1. We 
allow r(x) to be random as well, representing disorder 
in the hopping amplitudes, Jij- We could also allow J 
and u to be random, but this is less convenient and does 
not produce anything new. We write /i(x) = /i — e(x) 
and r(x) = ro + w{r) with disorder average [e(x)]au = 
[w(x)]a« — 0. For a given value of /x, the transition to 
superfluidity occurs with increasing J at a point Jc(/i). 
and we shall take S = J — Jdn)- 

Other representations of the problem serve equally 
well. The arguments we shall present are very general 



and do not depend on the precise form oi Cb- As seen 
below, the only common feature required is that time 
derivatives and chemical potential always appear in the 
combination [dr — fJ-)- In the coherent state functional 
integral formulation, for example, the Lagrangian has 
only a linear term, ip*{dT — /^)V'- The two formulations 
match up to irrelevant terms (in the renormalization 
group sense) so long as ^(x) does not vanish identically 
When /z(x) = the model ([T]) reduces to the well 
known random rod problem, corresponding to a classical 
(d + l )-di mensional XY-mode\ with columnar disorder 
m 3 This model maintains an exact particle- hole 



symmetry, and the Josephson scaling relation for k is 
therefore valid, but with z < d: see below. 



Twisted 0-houndary conditions [I4I are defined by the 
condition 



(2) 



where x = (t, x), e^, a — 0, 1, . . . , rf, are the space-time 
unit vectors, and is the dimension of the system along 
e„ (with Lo = /?)■ Thus '0(x) = e-*° ''~*"°'^V'(x), where 
^0 = 6*0//? and ko = {Oi/Li, . . .,6d/Ld), obeys periodic 
boundary conditions. We therefore expect ipo = (V'(x)) 
to be uniform and hence the order parameter (V'(x)) = 
giko-x-i-iujoT^^ to have a uniform phase twist. We now 
rewrite the Lagrangian in terms of ijj, obtaining 



Cb[iP; ^J.-i(^Q■,ro + J'kl]+'ko■P[^p]. (3) 



where P['(/'] = -iJ J d'^x dT[V'* Vt/'-VV^'*] is the total 
momentum. The simple shift, /i fi — iujQ, is guaranteed 
by the combination (dr—fi) for any model of dirty bosons. 

The free energy is defined by f^"'^" = {l3V)-^[lii{Z)]av, 
where the superscripts label the boundary condition. Us- 
ing ([3]), straightforward expansion of the free energy in 
powers of ko at = 0, 



ko 



(4) 



produces the usual momentum-momentum type correla- 
tion function expression for T. On the other hand, if 
ko = we have the obvious result /'^"(/i) = /(/i — it^o)- 
Analyticity of / within a given thermodynamic phase im- 
mediately implies the expansion 

Afr = f^"-f = luop + ^Ku;'^ + OiLol), (5) 

where p = ~df /dp and k = dp/dp- It is the comparison 
of equations ^ and ([5]) that motivates the identification 
of K with the temporal helicity modulus. Note, however, 
that the density p actually yields the leading term linear 
in ujQ. The reason this term appears in ^ but not in (j4|) 
is that Cb contains only quadratic spatial derivatives and 
is therefore even under space inversion, x — > — x. On the 
other hand P[ip] is odd under space inversion, implying 
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that /•'""o is an even function of ko. Moreover we see that 
kg breaks the inversion symmetry of In contrast, 

only if ^(x) =0 does £b possess inversion symmetry in 
time, T — > — T, which we call particle-hole symmetry 13]. 
This symmetry is already broken by ^, and the additional 
twist, loq, produces nothing new. These observations will 
be crucial to understanding the scaling of T and k. 

Let us consider first the Josephson scaling argument 
for the random rod problem, which then possesses both 
inversion and particle-hole symmetry. Classical intuition 
suggests in addition that k = in the disordered phase 
14| . As alluded to earlier, we then expect the lead- 
ing boundary condition dependence to enter via singular 
scaling combinations /cq? and Wq^t ■ Formally, this expec- 
tation is implemented via a finite-size scaling ansatz for 
the free energy density [1, , 



A/ 



ko^o 



where Aq and Bq are nonuniversal scale factors [15| . For 
convenience we take a hypercubical volume V = L'^, 
and the exponents are those of the classical random 
rod problem. The existence of a nonzero stiffness, i.e., 
a leading finite-size correction of order or /?~^, 

now requires that the scaling function ^q"'^" {x,y) « 
^d^oyZo'^o(^Q^x-^'^° +Qry^'^''°''°) for large x,y>0, yield- 
ing 



T 



(7) 



where 9^ = 9f + . . . + 9"^, implying the Josephson scaling 
relations vq = {d + zq — 2)i/o = 2 — ao — 2t'o and Vro — 
(d — ZQ)h'Q — 2 ~ aQ — 2zoi'o, and requiring in addition 
Qx,o Of 91q. 

Consider now the Bose glass to superfluid transi- 
tion. The density and compressibility of both phases are 
smooth nontrivial functions of ^, and a temporal twist 
now perturbs only slightly the particle-hole symmetry 
breaking term that is already in Cb- On the other hand, 
spatial twists still produce singular corrections, and at 
loq = we predict a finite-size scaling form [3] 



3l/zi/\ 



(8) 



with (x, y) k, R^x'^'^ 2)1/^21/ j^j, large y > 0, yielding 
T « ^('^-2)^B^'^(i?^/6l2)J^, v = (^d^z-2)v = 2-a- 
2^, and Rx oc 9^ as before. All exponents now refer to 
the dirty boson critical point. Because we have, as yet, 
imposed no temporal twist, there can be no 0(/?^^,/3^^) 
terms. 

Now, if a finite ujq is included, the only changes in ([8]) 
are that /i— icjQ everywhere, and boundary condition 
dependence of fa must be included. Thus 



A/ 



kflWo 



(9) 



where 5e = J — Jdl^ — i'^o) ~ ^ + i^oJ'c[^^)■ Most im- 
portantly, is the same function as that in ([8]), and 
therefore still yields no 0(/?~^,/3~^) terms. This implies 
immediately that the scaling function ^'^ft^o produces no 
direct contributions to p and k, which must therefore 
arise (a) indirectly from the dependence of Sg , and (b) 
directly from the analytic part of the free energy. The for- 
mer couples derivatives with respect to ^, or equivalently 
cjQ, to derivatives with respect to 5, producing the leading 
singular terms, Ks ~ I'^iT" [l3|- However, a = 2 — {d+z)u 
is very likely negative [3| , so this gives a vanishing contri- 
bution at the critical point. Therefore the main contri- 
bution is analytic in origin. Taking kg = 0, we may write 
/a(J,Ai) = -Pc{J)[lJ'- f^c{J)] - ^Kc{J)[fi- p.c{J)? + ■ ■ 
expanded for convenience about the transition line /J,c{J), 
and we immediately obtain from ^ a finite compress- 
ibility right through the transition, with the exponent z 
nevertheless undetermined 18, loj . 

There is a second approach to the derivation of the 
Josephson scaling relation, through scaling of the two 
point correlation function in the superfluid hydrody- 
namic regime, which we now discuss. Long wavelength, 
low frequency fluctuations are known to be governed by 
the effective Gaussian Lagrangian [1, 0] 

^^^ = \j d^'x J^^ dr[T\V4>\^ + K{dr4>)% (10) 

where (/)(x, r) is the coarse-grained phase, related to 
the coarse grained order parameter field via V-'(x, t) = 
"00 exp[i0(x, t)], where the bulk order parameter, V'o ~ 
\S\^ , near criticality. The small k, a; form of the Fourier 
transform of the two point correlator, G(x — x', t — r') = 
[{i;*{x,T)iP{x',T'))]ay, is then 



(11) 



Normally it is assumed that G obeys the scaling form 

G(k,w) w C\S\'''g{Dk\6\-'',Eu;\S\''''), (12) 

where 7 is the susceptibility exponent and C, D, E are 
again nonuniversal scale factors. If one naively matches 
([TT|) and (fT2|) one concludes that g{x,y) « (512^^ + 
g2y^)^^ for small x,y, and hence that 



(13) 



+ fa{J, ro + J^l, fJ- - iwo) - faiJ, ro,fi), 



Using the scaling relation a + 2/3-1-7 = 2 we obtain again 
the Josephson relations v — 2 — a — 2v, Vr =2 — a — 2zv. 

This argument, however, is just a disguised version of 
the free energy argument: £hd is obtained by assuming 
that the energetics describing global phase twists also 
describes slowly varying local phase twists. Thus, locally 
we replace ko by Vrf< and by dr(\>., then integrate over 
all of space-time [20[ . Since k arises from the nonscaling 
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part of the free energy, it is unlikely that it can now arise 
from the scaling part of the two-point function. Rather, 
we must carefully reconsider the scaling ansatz for G. 
Our proposal (for which we have no detailed theoretical 
support at this stage), is that 

|VolVG'(k,w) « H\6\'^'°'T{Dk\5\^'' ,Euj\5\-'"') 

+ ra(k,c^), (14) 

where is analytic. This self-energy scaling form is 
very similar to ([9]). Matching with pT|) . we again as- 
sume that T arises from the scaling part, k from the 
analytic part. Thus r(a;, y) w 71 cc^ -I- 722/^, for small x, y 
while Ta{k,Ljj) « koj"^ for small k,u!. The term yields 
T~|(5|'^. li z < d the term is subdominant to the 
analytic term, which yields n finite, as required, leaving 
z undetermined [21), l2^] . If we consider only static corre- 
lations, w = 0, standard scaling is recovered without any 
unusual analytic corrections. 

To summarize, the original scaling ansatz for the Bose 
glass-superfluid critical point, in which temporal twists 
produce a relevant symmetry breaking perturbation to 
the Lagrangian, and therefore scale in the combination 
wqCt-i is not supported theoretically, and may explain vi- 
olations oi z = dm. d = 2 Q. 
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